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Abstract. Particle flows in spatial networks are susceptible to congestion. In this 
paper, we analyze the phase transitions of these networks to a state of congested 
transport and the influence of both topology and spatial dynamics on its emergence. 
We systematically show that the value of the critical loading rate at which congestion 
emerges is affected by the addition of spatial dynamics, changing the nature of this 
transition from a continuous to a discontinuous one. Our numerical results are 
confirmed by introducing an analytical solvable framework. As a case of study, we 
explore the implications of our findings in the San Francisco road network where we 
can locate the roads that originate the congested phase. These roads are the spatially 
constrained, and not necessarily those with high betweenness as predicted by models 
without spatial dynamics. 
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1. Introduction 


Flow networks are inherently liable to congestion. The ability of these networks to 
handle demand at reasonable levels is crucial as otherwise a congested phase of transport 
affects the performance across the entire network. Therefore, it is of prime interest to 
analyze how and where networks begin to undergo a transition to a congested state 
and the dynamics of its response. Dissecting flow patterns is essential to address this 
problem. In this context, flow of data packets in the Internet is well understood, 
as analyses of their traffic dynamics and phase transitions are abundant [l]-j6]. The 
transition point to congestion is in this case well established through analytical solutions 
and simulations 00- To interpret the role of a network backbone for managing 
flows, optimal paths and minimum spanning trees have been studied 12 14 . The 


most relevant metric to determine the vulnerability of internet flows is the betweenness 
centrality for shortest paths, because it determines the critical element generating 


congestion in the network 15,16 


All these previous studies have overlooked the role of space and time in particle flows 
in the networks; a reasonable assumption in Internet applications. In this case the travel 
time between nodes is negligible and queues accumulate in the nodes which have limited 
capacity to process the packets and congestion occurs. In contrast, in transportation 
networks travel time of vehicles or individuals is crucial. Queues are formed on the links 
which have limited spatial capacity, and in turn their velocity depends on the density 
of travelers on each link. In this context, transportation research focuses on capturing 
traffic flow by making use of the fundamental diagrams that empirically relate flow, 
density and speed in road segments and utilizing them in macroscopic link models 


17,18. The cell transmission model 19,20 and the simple point queue models 21-23 


are well established among such traffic flow models. Alternatively, cellular automata 


models for vehicular traffic 24,25 have also been used to mimic traffic flow behavior 
along with many other discrete stochastic models 26,27 . In none of these cases, the 


interaction of the spatial dynamics with the network topology have been addressed. Our 
goal here is to build on the simple model of the Internet by adding temporal and spatial 
dimensions of particle flow and hence construct a framework to analyze congestion in 
spatial networks from a network science perspective. The proposed framework can be 
extended to flows in other kinds of spatial networks 28-30 , and more importantly 


could open new avenues of the research on urban road networks that would go beyond 


modeling and topological analysis for the statistical physics community 31-33 . 


2. Models 


We begin by recalling the scheme in BSIZH pi] and refer to it as the internet model 
(IM): the network is loaded with R identical particles at each timestep t with randomly 
assigned origins and destinations. A fixed shortest path routing table guides particles 
towards their destination. Nodes can transmit as many particles per timestep as their 
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outflow capacity , C, and travel between two nodes takes a unit timestep. Queues of 
particles form at the nodes, and they can grow infinitely large. Particles are exempt 
from joining the queue at their destination and are removed from the system upon 
arrival. The network response is measured by the order parameter H (9]: 

(A W) 


H (R) = lim 
v ' t^T RAt 


( 1 ) 


where W denotes the number of particles in the system, (AIT) is the average change 
in the number of particles still in the system after a timestep, At is the unit timestep 
and T is the length of the simulation. Figure [I] (a) and (b) depict the IM. For low 
values of R, the network reaches a rate of particle arrival equal to the loading rate. W 
remains constant and consequently H = 0. Conversely, if R exceeds a certain threshold 
R c , a linear increase in W with a slope of H is observed due to excessive queueing. This 
behavior maps a second-order phase transition to congestion. 

Spatial networks, being embedded in two-dimensional space, give rise to three 
heterogeneities that need to be captured. First is non-uniform travel times. This 
differs from the Internet where data packets hop from one node to the next in a single 
timestep. Second, these networks carry flows along the links. A particle in a spatial 
networks has a specific position on the link it is traveling on. The third source of 
heterogeneity is a consequence of this: particles occupy physical space and gradually fill 
the segment. Among the various models of traffic flow in the transportation literature 
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Figure 1 . (a) Number of particles W versus time t for the Internet model with 

different loading rates for the shown simple network with node and edge betweenness 
values mapped as node size and edge width, respectively, (b) A second order transition 
at the critical value R c = 7 for this network, (c) A schematic representation of the 
SPQM for a link with C = 2, r = 3 and V = 7. At t = 2, the link reaches its volume 
capacity and no incoming particle is allowed entering the link. 
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aimed to address these issues, the point-queue model (PQM) used in 21-23,34 is an 
adaptation of the IM that shifts flow from nodes to links, incorporating the non-uniform 
travel time distribution and thereby making the flow analysis very similar to that of the 
internet. Particles traverse the link freely by hopping through r slots of unit travel time 
to join a queue at the end of the link from which they will be discharged at the outflow 
capacity. The total travel time consists of the free travel time and the delay, namely, the 
timespan between the particle entering the queue and exiting it. The spatial point-queue 
model (SPQM) incorporates a single additional constraint to the PQM: every link has 
an upper limit for the number of particles it can hold at once. We will refer to this value 
as the volume capacity of the link, V. Links cannot accept any new particles when they 
reach their volume capacity, as illustrated in Figure [I] (c). This additional constraint 
has a crucial effect on the nature of the network response, as links at volume capacity 
clog upstream links and cause them to succumb to congestion as well. This spreading 
of congestion occurs at rates that depend on the loading rate and the network topology 
along with specific link properties. The rate of particles unable to travel determines the 
speed with which the congestion spills, which makes the spreading process non-binary 
unlike traditional spreading models in the literature. Although the SPQM share some 
aspects with several directed percolation models, the movement of non-identical particles 
along predetermined spatial shortest paths with non-binary spreading is uncommon and 


therefore relatively unstudied 35,36 


3. Results and discussion 

The critical loading rate has been shown |7j to be equal to Rff 1 = N(N — 
1) (Tmax/-Bmax) where N is the network size measured by the number of nodes, B ^ ax is the 
maximum node betweenness and C max is the outflow capacity of this node. This relation 
arises from the fact that inflow to a node is proportional to its betweenness centrality 
. At R CJ the inflow is equal to the outflow at the node with the minimum CAax/Anax 
value. For the PQM, we adjust this equation by replacing the node betweenness by 
an edge betweenness value, = N(N — l)(C' max /R^* x ), where £>^* x refers to the 

maximum modified betweenness. This modification is necessary since particles do not 
join the queues in the final links, which should be omitted from the edge betweenness 
calculations. 

Next, we introduce a framework to analytically calculate the entire transition curve 
to congestion. For R > R C1 particle inflow at certain elements will be larger than the 
outflow. We define R l c as the critical loading rate specific to element i. For large R , the 
outflow of congested links are maximized to capacity, which in consequence affects the 
inflow to the links downstream. To account for this we define the delay factor, Di(R), 
referring to the fraction of paths through i that are not suffering from delay as, 

D.(R ) = U(K ~R)+ j%;H(R - R% (2) 

where 'H(x) is the Heaviside step function and Ci and IfiR) are the outflow capacity 
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and the inflow of element i for loading rate R. Di(R) = 1 suggests no congestion for 
element i, whereas lower values indicate levels of congestion. Using this definition, the 
inflow Ii(R ) at a specific loading rate R can be quantified as, 

Sfcerp) II jek RDj(R) 


Ii(<W> N(N - 1 ) 

H(R) = - Ci)H(Ii(R) - Q), 

i£Af 


(3) 

(4) 


where T(z) is the set of paths passing through element i. Eq.Q accounts for all the 
delay factors of the elements upstream of element i by going through the shortest paths. 
Eq. (J2]) and (|3]) form a set of coupled equations that can be solved to obtain the inflows 
for every element. H(R) is obtained by summing all positive values of Ii(R) — Ci. 
Figure [I] (b) also reveals the exact solution for the simple network around the critical 
point with the dashed curve H = (7/42) (1 — R c /R). In order to test different network 
topologies and examine the effect of space on criticality, we use a non-periodic lattice 
as a substrate and rewire each edge (i,j) with probability p to a new destination j* 
chosen with probability proportional to d(i,j*) a where d(i,j) denotes the Euclidean 
. Figure [2] reveals that the simulations and analytical results perfectly 
coincide for transitions in both a two dimensional non-periodic lattice of iV = 1225 and 


distance 38,39 





loading rate R 

Figure 2. Transitions of a non-periodic lattice as estimated by simulations and 
by analytical solutions for (a) the IM and (c) the PQM. Transitions of a completely 
rewired network with p = 1 for (b) the IM and (d) the PQM for varying tendencies a 
to have long distance links. Analytical solutions are consistent with simulations. 
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C — 4, and its completely rewired instance. Critical rates decrease for smaller p as 
newly introduced shortcuts have higher betweenness values. As p increases, the value 
of a and its effect on network topology becomes more pronounced: rewired links in 
networks with lower a values are more localized and therefore maximum betweenness 
values are higher in these networks. Consequently the increase in H is sharper for lower 
a. Figure [2] (b) and (d) show sharper transitions and lower R c for decreasing a values 
for network instances with p = 1. 

At R Cl the link that triggers congestion, also referred to as the critical element, is 
expected to fluctuate between free flow and congested phases. Figure [3] (a) exhibits the 
frequency distribution of the timespans at which this most critical element operates at 
its outflow capacity, as an indicator of the temporality of the phase transition. Results 
show that these timespans follow a power law with exponents of —0.58 ±0.04 for the IM 
and —0.48 ± 0.04 for the PQM, independent of the network topology. SPQM exhibits a 
different behavior. For low volume capacities links tend to fill up, causing links upstream 
to fill as well. Fluctuations cause a gridlock , a condition where all elements of a cycle 
are completely filled and hence flows stop. In case of a gridlock, the order parameter 
increases very sharply. In Figure [3] (b), for varying volume capacities we measure the 
average number of timesteps it takes for a gridlock to occur, t 9 , normalized by the 
length of the simulation T. ft can be observed that for lower V, gridlocks are observed 
relatively quickly. For large volume capacities, PQM and SPQM have the same R Cl 
which we will refer to as the PQM-limit. In either model, at steady state, queues are 
not expected to be necessarily empty but rather steady in their size. If the volume 
capacity of a link is smaller than this steady state queue size, particles will be blocked 
in the upstream link which consequently may suffer from a decrease in its outflow due 
to this clogging effect. This suggests a lower critical point for a network with active 
volume capacity constraints, correspondingly f? p P M > RSPQM. Therefore, to realize the 
congestion-free transport to the fullest, the PQM-limit should be aimed. Figure [3] (c) 
illustrates the response of a non-periodic lattice with N = 625 for decreasing levels of 
volume capacity and reveals that both the critical loading rate and the nature of the 
transition is affected by V. Figure [3] (d) shows the effect of the volume levels on R c . In 
point queue models, the number of particles on a link is the sum of those in the queue, 
which was shown to be proportional to the modified betweenness, and those that are 
traveling. The expected inflow to a single slot of travel on a link is proportional to 
the actual betweenness of that link, hence the expected number of particles traveling 
is tB e . Therefore for small V and strictly deterministic inflows, a link reaches its 
maximum volume when V = R c (tB e ± /N(N — 1). The inflow to the link is 

proportional to the actual routing betweenness so the critical point can be expressed as, 




where B E is the edge betweenness. Consequently a linear increase in R c is observed up 
to the point where V = C(tB e + B E *) max /B E l x . However, the stochastic nature of the 
model, along with the first-order nature of its transition, causes fluctuations that force 
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Figure 3. (a) Distributions of the timespans through which the congested element 
operates at its outflow capacity at R c for the IM and PQM, following a power law. 
(b) The average number of timesteps normalized by T before a gridlock is observed, 
versus the volume capacity, (c) Change in the nature of the transition to congestion 
of a non-periodic lattice of size N = 625 with specified levels of volume capacities, (d) 
The critical loading rate R c for varying V/t values. R c settles to the PQM-limit as 
V/t increases. 


R c /C below this analytical bound as it converges to the PQ-limit. 

Figure [4] depicts the transitions for the PQM and the SPQM for the San Francisco 
road network with N = 1152 and an average degree of 3.2. The network is discretized 
by unit travel times of 10 seconds. Outflow capacity of a road segment is obtained by 
using the speed limit and the number of lanes. Volume capacities are estimated for every 
road segment assuming that the volume capacity is reached when speed drops to half of 
the speed limit. Under these assumptions, the PQM-limit is not reached as = 40 

(14400 vehicles/h) and /?y p Q M — 30 (10800 vehicles/h). To capture network response 
in the SPQM, segment volumes are recorded at different time steps of the simulation 
t = 720 (2h), t = 1440 (4h) and t = 2160 (6h) for R = 36, slightly above the critical load. 
Figure [4] illustrates the network response by mapping road segments that have reached 
80 percent of their volume capacities in the given time periods by color. Congestion 
originates from an artery leading to the downtown area and anisotropically spreads to 
other regions. After 6 hours, most network elements are suffering from congestion. 
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Figure 4. A congestion map of the San Francisco road network for R = 36 vehicles 
per timestep. Colors represent the times within which the road segments become 
congested as a consequence of the spillover. The black circle denotes the origin of 
congestion, (inset) PQM and SPQM transitions for the San Francisco road network, 
where R^ QM = 40 and Rf PQM = 30. 


4. Conclusion 


Our findings suggest that in the SPQM, the critical road segment cannot accept incoming 
vehicles as it saturates to its volume capacity before it reaches its outflow capacity. This 
outcome can be traced back to the inherence of congestion in downtown areas: cities 
with high population densities have concentrated spatial demand distributions, which 
result in the inadequacy of urban space to accommodate such concentrated flows. This 
work is a step further on a systems analysis applied to congestion in roads. In further 
studies, population and facility distributions can be modeled. In the applications domain 
an open question is to know how the studied transition to congestion is influenced by 
introducing empirical origin-destination matrices that represent the real population’s 
travel demand 40, 41 and how network topology contributes towards reducing- or 


maybe aggravating- congestion. While in the domain of phase transitions, analyzing 
the set of critical exponents using directed percolation as a benchmark case remains as 


an intriguing and open question 35,36 . 
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